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Abstract 

If P is a differential operator with constant coefficients, an identity is derived to calculate the 
action of e p on the product of two functions. In many-body theory, P describes the interaction 
Hamiltonian and the identity yields a hierarchy of Green functions. The identity is first derived for 
scalar fields and the standard hierarchy is recovered. Then the case of fermions is considered and the 
identity is used to calculate the generating function for the Green functions of an electron system in 
a time- dependent external potential. 

Keywords: 81V70 (Many-body theory), 16W30 (Coalgebras, bialgebras, Hopf algebras), 13N05 Deriva- 
tions 

1 Introduction 

In spite of immense progress in the calculation of the properties of materials, some of them are still beyond 
our computational power. Systems containing both localized and delocalized electrons often belong to 
these hard problems. A particularly conspicuous example is the optical spectrum of molecules and solids 
containing transition metals. Brute force calculations involving a large number of configurations do not 
give satisfactory results for the colour of such familiar materials as blood or grass In the solid 
state, non-perturbative Green function methods using the Bethe-Salpeter equation proved accurate for 
the colour of semiconductors [2112] E| but not for transition metal compounds. On the other hand, the 
very simple ligand field formalism can be used to calculate the position of the spectral lines jS], but not 
their intensities. The success of non-perturbative Green function methods for semiconductors and of 
ligand field theory for the spectral line position makes it desirable to unify the Green function and the 
ligand field theories. The present paper describes a tool developed for this unification. 

The main assumption of the Green function theory is that the eigenstate of the interacting system 
can be obtained as the adiabatic evolution of the ground state of the non-interacting system and that 
the latter can be written as a Slater determinant. A similar assumption is present in the Kohn-Sham 
approach to the density functional theory. On the contrary, the effectiveness of the ligand field approach 
comes from the fact that several initial states are taken into account. The ground state of the interacting 
system is then obtained by diagonalising an effective Hamiltonian representing the effect of the electron- 
electron interaction on these initial states 0. The starting point of our unification is therefore to use 
non-perturbative Green function methods to set up an effective Hamiltonian for a small number of initial 
states. The standard Green function method is recovered when the number of initial states is one, the 
ligand field method is obtained as the first approximation of a hierarchy of equations for the Green 
functions 

To develop this theory, we used a mathematical tool initiated by Schwinger [§| and that was taken up 
by a few authors |1(J1 1111 IT2] . This approach is known as the theory of Green functions for open shells 
|l,'i| and its development was rather slow because its combinatorial complexity is much larger than in the 
usual many-body theory. From a graphical point of view, the complexity comes from the fact that the 
Feynman diagrams must be drawn not only with the usual two-leg Feynman propagator but also with 
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2fc-leg "propagators" where k goes from 1 to the number of initial states. These many-leg propagators 
describe the correlation between initial states. It turns out that quantum group methods are well suited 
to tame this complexity and help deriving an identity that yields a hierarchy of Green functions for open 
shells |T3|. More precisely, if P is a differential operator with constant coefficients, the identity gives a 
closed formula for the action of e p on the product of two functions fg. 

For motivation, we show how this identity is used in many-body theory. The Green functions 
G(x\, . . . ,x n ) of an interacting system can be obtained from a generating function Z{j) by functional 
derivatives with respect to the external source j(x) (^H] p. 212) 



Therefore, in principle, Z(j) provides a complete information on the system. The generating function 
Z(j) is usually given by an equation of the type Z(j) — e p e w , where P is a polynomial in 5/5j(x) (the 
functional derivative with respect to the external source) and W° is a polynomial in j. In standard (i.e. 
closed shell) Green function theory, W° is bilinear in j: 



where G°(x,y) is the Green function of the non- interacting system. For the Green functions of open 
shells, the function W°(j) is a polynomial in j of degree larger than 2. This polynomial describes the 
initial state of the non- interacting system [lfi| . When we take the functional derivative of Z(j) with 
respect to j(x) and we denote SZ(j)/Sj(x) by Z'(J), we find 



z'U) = 




where Wq — SW°(j)/6j(x). The expression e p is an infinite sum of differential operators P n /nl acting on 
the product W^e . The role of the identity is to transform it into a finite sum of differential operators 
acting on e p e = Z(j). In other words, we obtain a non-trivial expression relating Z'(j) to a differential 
operator D acting on Z(j). This is indeed a hierarchy of Green functions, i.e. a relation between the 
one-point interacting Green function Z'(j) and some n-point interacting Green functions generated by 
DZ(j). Such a hierarchy is important because it provides non-perturbative equations that can be used 
to derive non-perturbative approximations |14| . In ref. |14|. the identity was given without proof and no 
detail was given to explain how the hierarchy follows from the identity. The purpose of the present paper 
is to fill this gap by proving the main identity and by deriving in detail some of its applications. 

The plan of the paper is as follows. After this introduction we show that the action of a differential 
operator P on a product of two functions fg can be described by a mathematical concept called a 
coproduct. Then we derive our main result which provides a way to calculate e p (fg). As a simple 
application of the basic identity, we derive the hierarchy of Green functions for standard quantum field 
theory of the scalar field. We also describe the changes that are required when we consider fermions 
instead of bosons. As an application, we derive in detail the generating function Z in the presence of an 
external time-dependent potential. We conclude with some possible extensions of the basic identity to 
non-commutative algebras. 

2 The coproduct of a differential operator 

In this section we introduce the coproduct of a differential operator with constant coefficients. We start 
with the case of a one-dimensional problem and generalize it to the d-dimensional case. But first, we 
define the tensor product of functions and of differential operators, that will be crucial for the definition 
of the coproduct. 

2.1 Tensor product 

If / and g are functions of a variable x, the tensor product of / and g, denoted by / ® g, is a function of 
two variable defined by (/ ® g)(x, y) = f{x)g(y). From this definition, it follows that the tensor product 



G[xi , . . . ,x n ) 



5 n Z 



Z Sj(xi) . . . Sj(x n ) ' 
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is linear: for any complex number A and for any functions /, g and h we have / ®(g + h) = f ® g + f ®h, 
(/ + g) ® h = f ® h + g ® h and A(/ ® g) = (A/) ® 5 = / ® (Ag). Moreover, we can multiply the two 
tensor products f®g and /'Ofir' by the rule (/ ® g)(f ® g') = (//') ® where (//')(&) = f(x)f'(x) 

and = g(y)g'(y)- 

If the variable a; runs over a d-dimensional space, we call differential operator a finite sum 

d m Qn d 

P- 2^ a ^-^dxY "'dxf' 

n 1 ,...,n d 1 a 

where a„ 1 ...„ £i are complex numbers. If P and Q are differential operators acting on functions of x, the 
tensor product of P and Q, denoted by P® Q is the operator acting on the tensor product of two functions 
I ®g by 

(P®Q)(f®g)(x,y) = ((Pf)®(Qg))(x,y) = (Pf)(x)(Qg)(y). 

With this definition, it can be checked that the tensor product of differential operators is linear and that 
we can define the multiplication of tensor products of operators by (P ® Q)(R ® S) = (PR) ® (QS). 
The tensor product is a familiar object in quantum physics to define an operator acting on a many-body 
state. For instance, if L is the angular momentum acting on a single particle, the angular momentum 
acting on two-particle wavefunctions is defined asZ/®l + l®L. 

2.2 The one-dimensional case 

If d = d/dx, the action of d n on the product of two functions of x is given by the Leibniz identity 

9 n (f9) = E f"V/)(3"- fe ff ). 



k=0 



In this definition, the functions / and g are not really important and look more like dummy arguments. 
The coproduct is a way to calculate d n (fg) without mentioning / and g. The coproduct is denoted by 
A and defined by 



Ad n = y~( r f)d k ®d n - 



For example 

Al = 1®1, 

Ad = <9®1 + 1®<9, 

Ad 2 = d 2 ®l + 2d®d+l®d 2 . 

More precisely, the relation between the coproduct and the derivative of fg is given by 

d n (.fg){x) = (Ad n )(f®g)(x,x). 

At this stage, we need to define more accurately the algebraic structure we are using. Let A be the 
algebra of differential operators with constant coefficients (i.e. the algebra of polynomials in the variable 
d). Thus P e A if and only if there is a finite number of complex numbers a n such that P = J2n>o a nd n - 
The product in A is the usual product of polynomials: if P = J2 n >o a nd n and Q — J2 m >o b m d m then 
PQ = ^2 n m>0 a n b m d n+rn . Thus, A is a unitary commutative algebra. The coproduct appears when the 
differential operators act on a product of two functions. It is a linear map A : A — > A ® A whose action 
on an element P € A can be written AP = ^2 t Ri <8> Si, so that P(fg) = J2i(Rif){Sig)- For example, 
for P = d 2 the coproduct is a sum of three terms described by R\ — d 2 , Si — 1, R2 — 2d, S2 — d and 
i?3 = 1, S3 = d 2 (note that, by linearity of the tensor product, we can also choose R2 = d, S2 = 2d). In 
general, we can use 



Hfg) = E a «E(fc)( 9fe /)( 9 "^), 

n>0 fe=0 ^ ' 
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to get 

AP = E a «E {f)d k ®d n - k . (1) 

n>0 fe=0 ^ ' 

In other words, the coproduct is just a way to denote which differential operators act on / and on g in 
P(fg). The introduction of such a concept may look rather pedantic at this stage, but it will allow us 
to obtain very general results, valid for partial derivatives and even for functional derivatives. Note that 
the fc-th derivative of d n with respect to d is nld n ~ k /(n — k)\ for n> k and for n < k, so that 

k>0 

where P^ is the fc-th derivative of P with respect to d. The concept of a coproduct is unfamiliar, but it 
is very useful in quantum theory. For instance, it was needed to calculate matrix elements of many-body 
operators and density correlations |16|. In this paper we shall see that it is also quite useful to derive 
non-perturbative equations in Green function theory. We hope that the example of the coproduct of a 
differential operator makes this concept understandable. For convenience, we shall replace the notation 
AP = J2t Ri ® Si by Sweedler's notation AP = Pil <S> Pi2, which is more common. 



2.3 The rf-dimensional case 

In the d-dimensional case, we put di = d/dxi and we define the coproduct of D = Y\t=i as 



fci=0 k d =0 

We can even define the coproduct of an infinite dimensional differential operator. If ri, . . . , r p are p 
points, the coproduct of the functional differential operator 

§m-\ \-n„ xrii fin p 

D <Jj(ri) ni ■ ■ ■ 6j(r p ) n * ~ Sjfay* ' ' ' 5j{r p ) n p 

is defined by 

jfeiH hfcp Jni— feiH \-Tip — hp 



AD = 



, _ , K kiJ " ' \k p J <5j(n) fc i . . .5j(r p ) k p 5j{vi) ni fcl . . . 5j(r p ) n p- k v ' 

K. \ — U r£ r, — U 



2.4 Algebra morphism 

Up to now, the coproduct is just a linear map from A to A. (8 A.. The coproduct becomes a powerful tool 
if it has a property called algebra morphism. The coproduct is an algebra morphism if Al = 1 <g> 1 and, 
for any P and Q in A, A(PQ) = (AP)(AQ) or, more explicitly 

A(PQ) = ^(PQ)il (g) (PQ)i2 = ^(PilQil) ® (Pi2Qi2) 

= (^Pil(g)Pi2)(^Qil(g>Qi2). 

To show that the coproduct of the algebra of differential operators is an algebra morphism, we consider 
for notational convenience the one-dimensional case. If we take P = d n and Q — d m , then PQ = d n+m 
and 



m+n / \ n 7n / \ / ' 

a(pq) = Ef n t m >^ 9n+m " = EEf!)f?)^ + ^^ 



k=0 1=0 



yn-\-m—k—l 



k=0 x 7 1=0 



^(l)d k ® d n - k ) ( J2 ( 7 ) dl ® dm ~ l ) = ( Ap )( A ^)' 



where we used the Vandermonde convolution formula X)fc+;=i (fe) (T) = (""j" 1 ) ■ The d-dimensional case 
is proved similarly. 
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2.5 Reduced coproduct 

For any element P of A, the reduced coproduct of P is defined by 

AP = AP - P <Z> 1 - 1 <gi P. 
For example: Ad = 0, A<9 2 = 2d ® 9, 

n— 1 / 

Ad™ = (fc) 5 *® 9 " - *' 

for n > 1. We denote the reduced coproduct of P by AP = P (1) ® P(2j- If the polynomial P is such 
that P(0) = 0, its coproduct AP contains the terms P ® 1 and 1 ® P. The purpose of the reduced 
coproduct is to eliminate these terms. As a result of this operation, the degree of P (i) and P (a) is always 
greater than and smaller than the degree of P. 



3 The main identity 

In this section, we prove the identity from which we can derive, among other things, the hierarchy of 
Green functions. Although this identity provides a powerful method to resum infinities of Feynman 
diagrams of the perturbation theory, its proof is very simple. 

Proposition 3.1. If A is a commutative algebra equipped with a coproduct A that is an algebra morphism 
and if P e A then 

A(c p ) = e ^ p (c p ®e p ), (3) 

where AP = AP - P®1-1®P and 

c ap = l0l + £_ ( AP r = i i + ]T_(£p a)0 P (2) ) 

7i—l n.— l 

Proof. By linearity of the coproduct, we have 

oo 1 OO _. OO ^ 

A ( eP ) - A (£^") = £^ A ( pn ) = £^( Ap )" = eAP < 

n— n— n— 

where we used the fact that A is an algebra morphism to write A(P") = (AP)™. In c AP , the operator 
M = AP is considered as a linear operator from A® A to A® A defined by M(Q®R) = ^(PilQ)®(Pi2R). 
We decompose this operator M as M — AP + P®l + l<g>P and we use the fact that the algebra A is 
commutative to write 

e M = e AP+P®l+l®P _ c AP e P®l e l®P 

The proof is completed by noting that 

OO OC 

e p®i = (p®i)» = y^ _p™®i = e p ®i. 

n=0 n=0 

Similarly, e^ p = 1 ® e p , so that e Pm e 1<g > p = (e p ® 1)(1 ® e p ) = e p ® e p . □ 

Note that A is not an algebra morphism so that (AP)™ ^ A(P"). Thus, we will use a special notation 
to write the terms of (AP)™. Namely, 

(AP)™ = ]TPil'"®Pi2'". 

With this notation, we can write the action of e p on the product of two functions as 

oo 

o P (fg) = (c p /)(e p 9 ) + ^l^(P 1 l'"c p / )(P 1 2'"e p ff ). (4) 

n— 1 
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4 Applications 



We consider now some applications of the identity Q of increasing complexity. The first application is 
a simple proof that e ad f(x) — f(x + a), the second is a hierarchy of Green functions for closed shells. 



4.1 The shift operator 

For this very simple application, assume that you do not know Taylor's theorem and that you want to 
evaluate e ad f(x) for an entire function /. To prove that e ad f(x) = f{x + a) it is enough to show that 
e ad x n = (x + a)"- The proof will be recursive. For n = 1, d°x — x, d x x = 1 and d k x — for k > 1. 
Thus, e ad x = (1 + ad)x — x + a. Now take n > 1 and P — ad. We have AP = -P(S>l + l<8>-P, so that 
AP = and the identity gives us e p (fg) = (e p f)(e p g). Therefore, taking f(x) — x n ~ x and g(x) = x 
we have e p (x n ) = (e p x n ~ 1 )(e p x). The recursion hypothesis tells us that e p x k = (x + a) k for all k < n, 
thus e p (x n ) = (x + a)™" 1 ^ + a) = {x + a) n . 



4.2 The standard hierarchy 

In the quantum theory of the scalar field, the system is described by the action A(ip) = A (f ) — V(<p), 
with 

A {tp) = i J dxd ti tp(x)d> 1 tp(x) - mV 2 (4 

and V(tp) — J dxF(ip(x)), where F(tp(x)) is a polynomial such that F(0) = F'(0) = (usually <^ 3 (ir)/3! 
or ( j c 4 (x)/4!). The generating function of the Green functions of this theory can be obtained from the 
following formula (see ^H] P- 445) 

Z(j) — e - iV (- iS i) e i f dxd yi( x ) G °( x 'y)j(y) ^ 

where Sj is a short notation for the functional derivative S/6j(x), and where G (x, y) is the Green function 
of the free scalar field. 

A hierarchy of Green functions is an exact relation between interacting Green functions. To obtain 
it, we first take the functional derivative of Z(j) with respect to j(x): 

62 (ti = c -iV(-iSi) 5 Q i f dzdz'i(z)G°(z,z'li(z') 

5j{x) 5j{x) 



<n-«*)( / dyG°(x,y)j(y)e^ d z dz' j (z)G (z, z ') j (z' ) y 



Now we use equation gj with / = / dyG°(x,y)j(y), g = e2 / dzdz'j{z)G°(z,z')j{z') and p = _ iV (_ i5 .y 
We first calculate (e p /) = e- iV( -~ iS ^ J dyG°(x,y)j(y). From the fact that F(0) = F'(0) = 0, we see 
that no terms of F has a degree smaller than 2. Therefore, V(—i5j) J dyG°(x,y)j(y) — because 
J dyG°(x,y)j(y) is of degree 1 in j(x). This gives us 

e -iv(- iSj ) f dyG\x,y)j{y) = [ dyG°(x,y)j(y). 



Now we calculate the reduced coproduct AV. If F(ip(x)) is a polynomial in ip(x), F(—iSj) is a polynomial 



in Sj. Thus 



Av - I ^HyAv) ®w-=I d yY, ( -ir 6 ^ ® F(k) (- iS m)- 

J fc>o m i(v) J fc>o K - 



Therefore, 

\k 



AV = 



I *yY, { n!r s m® F{k) (- i6 i(y))-v®i. 



k>l 
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The term n — 1 of the main identity Q becomes 

H) fe 



^ fe>i 



The fact that / = J dyG (x,y)j(y) is of degree 1 in j implies that only the term k = 1 may give a 
non-zero contribution. Its value is 

dy(-i)(8 j{y) f)(F'(-i8 m )e p g) = -i J dyG°(x,y)F'(-i6 m )Z(j). 

Moreover, all terms coming from (Ay)™ with n > 1 in the main identity Q give no contribution because 
they are of degree at least n in 8j. This gives us the simple result 

' )/{J] = i [dyG°(x,y)j(y)Z(j)-i [dyG°(x,y)F'(-i8 m )Z(j). (5) 



5j(x) 

If we multiply this equation by (□ + m 2 ) we get the standard hierarchy (see ^H] P- 447) 

which is usually obtained by path-integral methods. However, note that equation JJJJ does not imply 
equation JJjJ but 

= iJdyG (x,y)j(y)Z(j)-iJdyG (x,y)F'(-i8 j{y) )Z(j) + 4>(x), 

where <j>{x) is a solution of the scalar wave equation (□ + m 2 )0(x) = 0. Thus, equation JSJ is a result 
stronger than JSJ. 

5 The fermionic case 

For applications to molecular or solid-state physics, we need to modify the equations to treat the case of 
fermions. The main difference with the previous case is in the definition of the tensor product. 

5.1 Tensor product 

The basic variables of the fermionic theories are the field operators ip( x ) an d i/j'(x) and the fermionic 
sources r](x) and fj(x). These variables are assumed to anticommute: for example rj(x)ip(y) — —ip(y)rj(x). 

The functions / will be polynomials in these basic variables. If / is a monomial in the basic variables, 
we denote by deg(/) the degree of /. For example, if / is a basic variable, deg(/) = 1. For / = ^t(r)^>(r), 
deg(/) = 2; for / = J f}(x)ip(x)dx, deg(f) = 2. The parity of a function /, denoted by |/|, is if deg(/) 
is even and is 1 if deg(/) is odd. Moreover, a function / is said to be even if |/| = and odd if |/| = 1. 
From the anticommutation of the basic variables, one can derive the following commutation relation of 
two functions / and g: gf — (— l)'^" 9 '/g. Therefore, an even function commutes with all functions and 
two odd functions anticommute. 

The differential operators are products of functional derivatives with respect to 77 or fj. The basic 
relations are 

5 S 

-r](y) = S(x-y), ——fj(y) = 0, 



8rj(x) ' 5rj{x) 

8 S 

-My) = S(x-y), ——-ri(y)=Q. 



5r)(x) ' 8rj{x) 

If / and g are functions in the basic variable with a definite parity, then the functional derivative of the 
product fg is given by the modified Leibniz formula 

8r](x) ^ 8rj(x)® ^ ^ ^ ^ 57](xY ^ 
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and the same equation for a functional derivative with respect to fj{x). Equation (JJJ is known as Leibniz' 
rule. 

The sources rj and fj anticommute, so the functional derivatives anticommute: 

8 2 S 2 



8rj(x)8r)(y) 8r}{y)8rj{x) 

The degree of a monomial P in 5/5r)(x) and S/Sfj(x) is denoted by deg(P) and the parity of P, denoted 
by \P\ is if deg(P) is even and 1 if deg(P) is odd. The action of a differential operator on a product of 
two functions is given by 

P(fg) = E(- 1 ) l/l|Pl2| ( Pll /)( Pl2 5) ; 

where the coproduct of terms of degree 1 is 

. 8 8 , 1 8 

A——- = - ; ; ® 1 + 1 < 



Srj(x) 8rj{x) 8rj(x) 

8 8 8 

A — ; = — ; 1 + 1'' 



8rj(x) 8rj{x) 8rj(x) ' 

and the coproduct of a product of two monomials P and Q is calculated from the coproduct of P and Q 
by using the fact that the coproduct is an algebra morphism: 

A(PQ) = (AP)(AQ) =^(-l)l Pl2 IIQ ll l(PilQil)® (Pi2Qi2), (8) 
where we have used the fact that the multiplication of P ® Q by R ® S is (-l^H^PP) <g> (QS). 

5.2 The main identity 

With the conventions described in the last paragraph, the main identity J3J is valid without change for 
fermionic sources if P is even (i.e. if |P| = 0). The action on a product of two functions is now 

oo 

e P (f 9 ) = (e p /)(e^) + ^i^(-l)l/H^ m l(P 1 l'"(c p /))(P 1 2'"(e^)). (9) 



n 

n=l 



6 Application of the fermionic formula 

The main application of the fermionic formula was the derivation of the hierarchy of Green functions for 
open shells ^31- It would be too long to give a detailed derivation of this hierarchy. So we consider the 
calculation of the generating function Z of the Green functions for a system in a time-dependent external 
potential. This generating function is useful in time-dependent density functional theory. We first describe 
some properties of fermionic sources, then we give the step by step derivation of the dependence of Z on 
•q and fj, using the main identity. Finally, we use this identity again to obtain the exact expression of Z. 

6.1 Fermion sources 

The S-matrix of a fermionic system with external sources rj and fj and action ^4 mt is 
S(fj, rf) = Texp ( - iA iut + i / fj(x)ip(x)dx + i / V 1 '{x)rj(x)dx) . 



For a nonrelativistic fermion, i/j(x) and ip (x) have a spin index. Thus, the sources are also two-component 
vectors and 

2 2 

fj(x)ip(x) = y]fj s (x)ip s (x), ip^(x)ri(x) = y]i/}l(x)r} a (x). 



To see how functional derivatives act with respect to the time-ordering operator, we first notice that 
the sources can be taken out of the time-ordering operator. For example, if x° > y 

T(fj(x)ip(x)ipi (y)r)(y)) = f){x)^(x)^ {y)r]{y) = fj(x)rj(y)tl)(x)^ (y) 
= v(x)v(y)T(^{x)^(y)), 

if a; < y° 

T(fj{x)ip{x)ip*(y)r)(y)) = ^ (y)i](y)fj(x)ip(x) = -fj(x)r)(y)ij} 1 '(y)i>(x) 

Thus, the functional derivative with respect to rj(x) or fj(x) commutes with the time-ordering operator. 
In particular |15| 

8S(fj,rj) t _ :rrLl js- iA ™ 



Srj(x) 



k=„=o = iT(i>{x)e- lA ), 



where the minus sign in the last equation comes from the fact that the functional derivative must jump 
over %Ip{x) to reach r)(x) in the definition of S(fj,ii). 

A standard result of the functional derivative approach |15II17| is that the interacting S-matrix S(fj, rj) 
can be obtained from the non-interacting S-matrix 

S°(fj,?]) = Texp(i J fj(x)ijj(x) + ^ (x)i](x)dx^j 

by the equation 

/oo ■ c ■ c 

where H lnt (ijj(r, t), tp^r, t)) is the interaction Hamiltonian. 

If |$o) is a non-degenerate eigenstate of Hq, the Green functions of the interacting system can be 
obtained from the generating function 

Z%r{) = (* |5(77,7,)|*o> = exp(-i / H^(—^, )dt)Z°(fj, rj), 



1 5r)(x) ' 5r)(x) ' 

where Z° = ($o|<S' | ( I > o)- If \®o) can be written as a Slater determinant, |$o) = b\ .. . b\ 0), where b\ is 
the creation operator of the one-electron state u n (x) of the free Hamiltonian, we have j!6j 



Z°(fj, 77) = (^ol^exp (i I r){x)ip{x) + ^(x)r)(x)dx) |* ) = e"^ n[x)G fovMv)****, 

where 

G°(x,y) = -i($ \T(l>(x)tf(j,))\$o) = ^(0\T(^(x)^( y ))\0)+ l Y / ^X k (y)- 



N 



k=l 



As a useful example, we consider the interaction with an external time-dependent potential v(x) so 
that A lnt = J dxip> (x)v(x)x/j(x). The generating function becomes Z = e p e w , where 

f S 2 
P = —i dxv(x)- 



8rj(x)5r)(x) ' 
W a = -i I f](x)G°(x,y)r 1 (y)dxdy. 



Our purpose is now to calculate Z. 



6.2 Dependence of Z on r] and 77 

For later convenience, we consider Z(X) = e p e xw ° . We first calculate the derivative of Z(X) with respect 
to fj(x). 



We can apply identity © with / = f G°(x,y)rj(y)dy, so that |/| = 1 and g = e xw ° . We first note that 
Pf = because / does not contain any factor fj. Thus e p f = f and since e p g = Z(X), equation (0 
becomes 



XfZ(X) -iXj^^ E(- 1 )' Fl2 " J| ( Pll '"Z) (^i2" l Z(A)) . 



<JZ(A) 

v 7 n— 1 



To proceed, we must calculate the reduced coproduct of P. We first have the coproduct 

f S S f A A 

AP = P <E>1 + I <Z> P - i / dzv(z) ® + i / dzv(z)- 



Srj(z) Srj(z) J 5r)(z) 8rj(z) 



5 S 

AP = -i I dzv(z) rr ® 7=TT +^ / dzv(z) " ® " (10) 



so that 

<5?7(z) <5?7(z) 7 <5r/(z) #77(2) 

The term / is of degree 1 in rj and in fj. Thus, the factor Pil' n f is non zero only if n = 1 and 
Pil' = 6/6r)(z). We find 

^ (A) - -iXfZ{X)-iX{-i) [dzv(z)(-l)- 6f 5Z{X) 



Sr)(x) ' J 5i](z) 8rj{z) 

SZ(X) 



iXfZ(X) + X / dyv(y)G°(x,y) 



6r)(y) 



If we write Z(X) = e w< -^ we get 



with the solution 



= -ix[d y G*(x,yM y ) + x[d y G^ X ,y)v(y) 6 -^., 
drj(x) J J drj(y) 



M1 ' Ai = -iA / dyG(x,y) V (y), (II) 



<5t?(:e) 

where G(x,y) satisfies the Dyson equation 



G(aj,y) = G°(x,y) + X dzG°(x,z)v(z)G(z,y). (12) 



In other words, G(x,y) is the Green function of the Schrodinger equation obtained by adding Xv to 
the free Hamiltonian. From equation (|llfl . we see that the function W(X) can be written W(X) = 
— iX J dxdyfj(x)G(x,y)r](y) + g(ri,X). The same calculation for 5Z(X)/Sr](y) shows that SW(X)/5r](y) = 
iX J dxfj{x)G{x 1 y). Thus, g does not depend on rj and 



W(X) = -iX I dxdyfj(x)G(x,y)r)(y) + g(X), (13) 
where g(X) is independent of fj and rj. 
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6.3 Calculation of 17(A) 

To calculate 17(A), we take the derivative of Z(X) with respect to A. 



OZ(X) _ e P^ W (^e-° 



dX V d\ 

We can apply identity © with / = W° , so that |/| = and g = e xw ° . We now have 

PW° = -i / dxvix) — — — -r = - / v(x)G°(x,x)dx = -tr(«G°). 
J dr][x)dr]{x) J 

Thus e p f = f — tr(i>G°) and since e p g = Z(A), equation becomes 

= W Z(X)-tr(vG )Z(X) + f^^J2( Pll ' nw0 )( Pl2 ' nZ W)- 

n—l 

To calculate the term n = 1, we use the reduced coproduct of P given in equation (|10|l and we obtain 

W" *Z(A) . f , , . W° <JZ(A) 



5rj(z) Srj(z) 



dxdyfj(x)G°(x,y)v{y)^^^ + J dxdyv(x)G (x 1 y)r 1 (y)^^. 
For the term n = 2 we first evaluate, using algebra morphism ijSJl 

(AP) 2 = C - i ( dxv(x) S <g> ~4 + i [ dxv(x) * 



£77(2;) ^(x) 7 <5?7(x) ^(x) 

( - i / dyw(y)—^— <g> -J-^ + i / dyw^)-^- <g> -j-^r) 
V 7 0T7(2/J drj(y) J drj(y) drj(y)/ 

S 2 

5i](x)Sr/(y) Sfj(x)6fj(y) 
dxdyv(x)v(y) - 

dxdyv(x)v(y) - 

dxdyv{x)v{y) - 



5r}(x)6r)(y) 


5r)(x)6rj(y) 


S 2 


S 2 


6fj(x)5i](y) 


6r)(x)5fj(y) 


S 2 


S 2 


Sfj(x)Sfj(y) 


6r)(x)5r)(y) 



W° contains one 77 and one fj. Thus, the only terms of (AP) 2 that give non-zero contributions are the 
second and the third, because the left hand side of the tensor product contains functional derivatives with 
respect to rj and fj. If we interchange the variables x and y, we see that these two terms are identical, 
thus the term n — 2 becomes 



i^(Pil'V°)(Pi2 ,2 Z(A)) = - j dxdyv{x)v{y) 



S 2 W° S 2 Z(X) 
Sr/(y)Sfj(x) Sfj(y)Srj(x) 

= i I dxdyv(x)v(y)G° (x, y)- 



8 2 Z{X) 
' 6fj(y)Srj(x) 



There is no term for n > 2 because W° is of degree 2. 
This gives us the equation for Z(X) 

,>Z{X) W°Z(X)-tr(vG a )Z(X)+ [ dxdyfj{x)G°(x,y)v(y)- Z ^ 



dX J Srj(y) 

dxdyv(x)G°(x 1 y)r](y) SZ ^ +i [ dxdyv(x)v(y)G° (x, y) ^ . 

drj(x) J drj{y)dr]{x) 
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If we write Z(X) = e w ( x > we obtain the equation for W(X) 

= -iJfj(x)G (x,y)r,(y)dxdy-tr(vG )+Jdxdyfj(x)G°(x,y)v(y) 6 -^ 

+ f dxdyv(x)G°(x,y)r](y)^^-+i f dxdyv{x)v(y)G°(x,y) SW ^ SW ^ 



8rj(x) J Srj(y) 8rj(x) 



+ iJdxdyv(x)v(y)G°(x,y)^^ 



The general form of W(X) is given by equation If we introduce it into the last equation we obtain 

-irjG°ri - tr(vG°) - iXf]G G vGr] 
-iX 2 fjGvG°vGr] - Xtr(vG°vG). 



BC 

-ifjGr>-iXfj—r) + g'(X) = -^G ^ - tr(vG°) - iXfjG°vGri - iXfjGvG ^ 
BX 



This gives us two independent equations: 

BC 

G + X— = G + XG°vG + XGvG + X 2 GvG°vG, 
BX 

g'(X) = -ti(vG°) - Xtr(vG°vG). 
Equation (|12|> enables us to simplify this into 

BC 

G + X— = G + XGvG, 
BX 

g'(X) = -tv(vG). 
The first equation is an identity because the Dyson equation l|12l) yields 

^ = G°vG + XG°v^-, 

BX BX 

so that 

It = GvG - 

BX 

The last equation to solve is simply g'(X) = — tr(uG) = —i J v(x)p(x; X)dx, where p(x; X) = —iG{x, x) is 
the charge density in the presence of the potential Xv. The solution of this equation would be g(X) —g(0) = 
g'(fj,)dfj,. If A = 0, Z(0) = c p l = 1. Thus, W(0) = 0. Moreover, for A = 0, G = G° and W = g(0). 
Thus, g(0) = 0. Now we prove that 

g(X) = -tr(log(l + A V G)). (14) 
By definition, the right hand side of this equation is 

-tr(log(l + AuG)) = ^LjV_tr((«G)"). 

n=l U 

If we take the derivative with respect to A we find 

_ at ,(log(l + A, G )) _ ^ ( „ rll[M1 ^ H) , t[( « ( , r , 

ri—l n— 1 

OO CO 

= - ^(-A)"tr((«G)" +1 ) + ^(-A)»tr((«G)" +1 ) = -tr(vG). 

n— n— 1 

Therefore, g(X) and — tr(log(l + XvG)) satisfy the same first order differential equation. Moreover, the 
boundary conditions are the same because — tr(log(l)) = = g(0). 
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6.4 Alternative expressions 

We can give alternative expressions for <?(A). Expanding G = G°+XG°vG we obtain G = G° J2^Lo(^ v G°) n 
so that 

oo 

g'(X) = -tr(«G) = -^A"tr(( V G°r +1 ). 

n=0 

Hence, 

r\ 00 vn+l 

5 (A) = / 5 '( M ) d/ i = -^— -tr(( V G r +1 )-trlog(l-A V G ). 
Jo t'o n+1 

A similar result was obtained by Sham in the framework of the density functional theory. The relation 
with equation 12} is obtained from G = G° + \GvG°, so that vG° = vG(l - XvG°). Therefore, formally, 

trlog(l - AvG°) = trlog(uG°) - trlog(wG) = trlog(wG°) - trlog(uG°(l + XvG)) 
= -trlog(l + XvG). 

We can manipulate these operators as if they were scalars because vG is a function of vG°, so that every 
term is a series in the single variable vG°. 

To give a last form of <?(A), we define the Euler operator A = J dxv(x) Sv S ^ x j , that counts the number of 

times v is present in an expression. In other words, it is easy to check that Atr((vG )") ~ rttr((vG ) n ) . 
More generally, A k tr((vG°) n ) — n fc tr((uG )") . Thus, we can write the somewhat formal identity 

00 ,„ +1 00 
,g(A) = -J2 ^tr(( V G )" +1 ) - -A- 1 Y, X n+1 tr((vG°) n+1 ) = -XA-Hr(vG). 

n=0 U n=0 

7 Conclusion 

The hierarchy of Green functions is one of the building blocks of an extension of the Bethe-Salpeter 
equation to open shells. The other main ingredient is an expression of the Green functions in terms of 
2-particle-irreducible generating functions. This expression is well known for closed shells but not 
for open shells. In the latter case, the Dyson equation relating the two-point Green function and the 
self-energy (i.e. one-particle irreducible Green function) is not valid. This is a bad news because the 
Dyson equation is required to use the Legendre transformation relating the connected Green function 
and the one- or two-particle irreducible ones. By a diagrammatic tour de force, Hall was able to obtain 
the equation corresponding to the Dyson equation for open shells [111 . It looks very complex but, by 
doubling the size of the Green functions, it is possible to give it a form similar to the traditional Dyson 
equation. This gives some hope that the Bcthc-Salpeter equation for open shells can be obtained, opening 
the way to an effective unification of the ligand field and Green function methods. 

In some applications, the operator P in e p belongs to a noncommutative algebra [201 - This is the case 
for example when P is a differential operator with non-constant coefficients. It is possible to generalize 
our main identity to noncommutative algebras. The main problem is to transform e AP into an operator 
e A acting on 1 ® e p . In other words, we have to solve the equation e A — e Ap e _1(8P . This can be done 
by using the Baker-Campbell-Hausdorff formula [5T] . 
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